Weakly-coupled Hubbard chains at half-filling and Confinement 
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We study two (very) weakly-coupled Hubbard chains in the half-filled case, and especially the 
situation where the intrachain Mott scale m is much larger than the (bare) single-electron inter- 
chain hopping t±. First, we find that the divergence of the intrachain Umklapp channel at the 
Mott transition results in the complete vanishing of the single-electron interchain hopping: This 
is significant of a strong confinement of coherence along the chains. Excitations are usual charge 
fcrmionic solitons and spinon-(anti)spinon pairs of the Heisenberg chain. Then, we show rigorously 
how the tunneling of spinon-(anti)spinon pairs produces an antiferromagnetic interchain exchange of 
the order of J± = t± 2 jm. In the "confined" phase and in the far Infra Red, the system behaves as a 
pure spin ladder. The final result is an insulating ground state with spin-gapped excitations exactly 
as in the opposite "delocalized" limit (i.e. for rather large interchain hoppings) where the two-leg 
ladder is in the well-known insulating D-Mott phase. Unlike for materials with an infinite number of 
coupled chains (Bechgaard salts), the confinement /deconfinement transition at absolute zero is here 
a simple crossover: no metallic phase is found in undoped two-leg ladders. This statement might be 
generalized for N-leg ladders with N=3,4... (but not too large). 

PACS numbers: 71.10.Pm, 74.20.Mn 



I. INTRODUCTION 

One-dimensional (ID) electron systems have attracted 
a great attention over the last years. The Hubbard chain, 
that is a nice prototype to describe ID conductors, sim- 
ply reduces to the so-called "Luttinger liquid" at suffi- 
ciently low energy 0] . Due to the restricted motion along 
one direction in space, perturbations easily propagate co- 
herently. This induces that spin and charge degrees of 
freedom get independent, and that all the low-energy ex- 
citations are collective modes - namely, longwavelength 
fluctuations of the charge- or spin density. No Landau 
quasiparticle type elementary excitations exist . 

Two-coupled Hubbard chains have been studied as 
a basic model which includes intrachain interaction U, 
longitudinal t and transverse t± hoppings 0. For re- 
cent and general reviews on this subject, consult Refs. 
§-0]. In this paper, we study the case of a weakly cou- 
pled Hubbard ladder at (and close to) half-filling with 
t± <C U <C t. It should be noted that this interesting 
situation has been rarely investigated in the litterature. 

The presence of intrachain Umklapps provides two 
competing energy scales || , the Zeeman-like band split- 
ting energy tending to delocalize (deconfine) particles in 
the transverse direction [01 



A cx t 1 / 1 '* 1 with T) cx U 2 , 
and the well-known single-chain Mott scale 
to oc exp — TTVp/U, 



(1) 



(2) 



that may rather induce confinement along the chains. In- 
fluenced by the plethora of novel phenomena predicted 



(postulated) at the connnement/deconfinement transi- 
tion in systems with an infinite number of weakly cou- 
pled chains (Bechgaard salts) {Hi] , it is naturally worth to 
consider such issues on a two-leg ladder where tractable 
calculations are indeed possible. 

When the splitting energy dominates - i.e. when the 
hopping amplitude increases much faster than the Umk- 
lapp channel(s) upon renormalization - that would corre- 
spond to the deconfinement phase. The result is a two- 
band model with four Fermi points. Rewriting the bare 
interactions - Do consult Eq.(||) - in the so-called band 
basis, we immediately recover the model that has been 
previously studied by Lin, Balents and Fisher in the limit 
of very large interchain hoppings A> (7 ]Tl| ] . The renor- 
malization group transformation scales the system to a 
special strong-coupling Hamiltonian with enormous sym- 
metry - the SO (8) Gross-Neveu model, that is integrable. 
This allows to conclude that as soon as A > m the 
two-leg ladder system is in the so-called D-Mott insulat- 
ing phase with spin-gapped excitations and "preformed" 
short-range d-wave pairing correlations |T^ | . 

The splitting of the two bands affects spin-charge sep- 
aration: Cooper pairs and magnons appear as natu- 
ral excitations both at the two-band Mott scale M w 
t± exp — ttvf/8U that is of the same order as m. Dop- 
ing such a spin liquid liberates "hole-pairs" with d-wave 
pairing, restoring partially spin-charge separation 
There is an extended massive 7r-mode with S0(6) sym- 
metry, containing remnant charge and spin degrees of 
freedom [p^| . 

Then, we study this model of interacting electrons hop- 
ping on a two- leg ladder, focusing on the behavior at 
half-filling mainly in the opposite limit where the intra- 
chain Umklapp vertex is flowing first to strong couplings 
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- i.e. for very small interchain hoppings. First, we show 
that the divergence of the Umklapp channel at the Mott 
transition makes the single-electron interchain transfer 
completely vanishing, resulting in perfect confinement of 
coherence along the chains. We insist on the fact that 
the Kosterlitz-Thouless like transition is then associated 
with a jump both in the Luttinger liquid exponent and 
in the single-electron interchain transfer. 

We may explain the entire suppression of the latter at 
the Mott scale, as a consequence of the strong breakup 
of physical electrons. The system can be also seen as a 
superposition of two degenerate bands. A noticeable fact 
is that in absence of Umklapps i.e. away from half-filling, 
this condition of confinement is never satisfied for on-site 
(weak) interactions JlJ]. However, this could happen at 
zero temperature in presence of long-range interactions 
along the chains Jf5|,[l6| (or still low-frequency Holstein 
phonons), or possibly in very restrictive models extended 
to an infinite number of coupled Hubbard chains fl7f . 

At the opening of the Mott gap, low-lying excitations 
located mainly on the chains are usual empty/doubly oc- 
cupied sites and spinon-(anti)spinon pairs. Then, as a 
result of spin-charge separation arising in the confined 
phase, we show rigorously how the tunneling of spinon- 
(anti)spinon pairs give rise to an interchain Heisenbcrg 
exchange coupling of the order of J± = t± 2 /m > [Q. 
This has been previously neglected in Ref. [Q. The 
result is equivalent to two weakly coupled Heisenberg 
chains p(J. The ground state is a disordered spin liq- 
uid with a prominent Akp Charge Density Wave (CDW), 
like in the deconfined limit. It is worth noting that the 
spin gap w J± rejoins the charge gap in the entrance of 
the D-Mott phase. 

To summarize: The confinement /deconfinement tran- 
sition arising in undoped two-leg ladder systems is a pure 
crossover UW. No metallic phase is found. This state- 
ment might be generalized for a N-leg ladder with finite 
N (N = 3,4...). Conclusions are really different than for 
an infinite number of coupled chains where the decon- 
finement transition |^2|,|| - i.e. when formally A ~ m - 
would be a priori responsible for the difference of behav- 
ior between the TMTTF (ID insulator) and the TMTSF 
(strange metal) Bechgaard salts jl0| . 

The plan of the paper is as follows. In Sec. II, the 
model and formalism are introduced. In Sec. Ill, based 
on Renormalization Group arguments and bosonization 
tools, we define the confinement phenomenon. In Sec. 
IV, we study in details the spinon pair hopping: The 
system behaves as a pure spin ladder. In Sec. V, doping 
effect on the confined state and incoherence of t± away 
from half-filling are carefully studied. In Sec. VI, we for- 
mulate a brief summary of our results and make compar- 
isons with Bechgaard salts. In Appendix A, (Abelian) 
definitions of spinon- and charge excitations below the 
Mott scale are given. In Appendix B, "anyonic" charge 
excitations of Luttinger liquids are reminded. 



II. MODEL 

We start with the following model. The kinetic energy 
takes the form TLkin =H + T~t± 



t± y2 <T 2a (x)d la (x) + H.i 



(3) 



The indices j = 1, 2 denote here the chains, and a =|, J,. 
Focusing on electronic states near the Fermi points, one 
can expand: 



d ja (x) = d +ja {x)e 



ikpx 



d-i 



[x e 



-ikpx 



(4) 



At half-filling, one must equate kp — ir/2. We put the 
bare short-distance cutoff, a — I. Then, the Hubbard 
four-Fermion interaction can be expressed in terms of 
currents, defined as [jll 12 



(5) 



a, a 



/ dpj a € aa ' dpj a ' , 



and in the following p = ± denote respectively right and 
left excitations, a denote Pauli matrices and e Qa < is an- 
tisymmetric: e aa > — —e a / a and e-j-j, = I. Precisely, the 
set of marginal momentum (and non-) conserving four- 
Fermion interaction reads: 



Hint — 9cJ+jjJ- 



33 



9 s J 



+33^-33 



+33^-33- 



(0) 



Note that g c and g s describe charge- and spin backscat- 
tcrings respectively, and g u the intrachain umklapp pro- 
cesses. The bare interactions are of the order of U. 

Now, it is appropriate to use an Abelian bosonized 
form for fermion operators. These transform as [§j3| 

d pja = exp [p{<S>jc + a®js) - (% c + a9 js )} ). (7) 

a = ± for spin up and spin down, respectively. Re- 
markably, electron spectrum yields spin-charge separa- 
tion. Absorbing the interaction g Cl that is not affected 
by a rescaling of the short-distance cutoff, in the charge 
part of H results in the so-called Luttinger model 



H 03 = ^j dX ^ 



(8) 



dx&jc = {pjc~Po) measures fluctuations of charge density 
in each chain, and VQj C is the conjugate momentum to 
<j>j C . All the interaction effects are now hidden in the 
parameters u (the velocity of charge excitations, uK = 
vp with the Fermi velocity vp = It sin hp) and K (the 
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Luttinger liquid (LL) exponent controlling the decay of 
correlation functions, K = 1 — g c /i:vF)- K < 1 means 
repulsive interactions. 

The free spin Hamiltonian yields the same form replac- 
ing $j c by and the Luttinger exponent K by K s = 1 
due to the requirement of SU(2) invariancc. 



III. CONFINEMENT DUE TO UMKLAPPS 

It is first suitable to write the main Renormalization 
Group (RG) equations. 



A. RG-flow analysis 



By simple scaling arguments, one gets: 
dlnt ± 3 1,„ 1, 

-^r = 2-A {K+ K ) - 



(9) 



Anisotropy between space and time induced by Umklapp 
scattering results in a strong renormalization of the LL 
exponent and the charge velocity ]23[|: 



^ = -C x {g u K)\ 

du n 2 V 

-77 = -C 2 g u uK. 
dl 



(10) 



For a complete derivation of such equations, do consult 
Rcf. p4| ]. The / describes renormalization of the short- 
distance cutoff a(l) = expl. At half-filling, usual Bessel 
functions Jo and Ji are non-oscillating and have been 
included in the positive constants C\ and C2 p3j- If 
one is at finite temperature T, then the renormalization 
procedure must be stopped at lengths a{l) comparable to 
the thermal length u/T, that means I ~ ln(l/T). 

Starting with a sufficiently small (bare) interchain hop- 
ping, then the Umklapp channel g u is flowing first to 
strong couplings at a critical scale T c that is the single- 
chain Mott scale m given by Eq.(||). For one- loop RG 
equations, indeed one gets a critical driving parameter 
l c = avp/U where a is of the order of unity. Although 
the singularity at finite l c is an artifact of the one-loop 
calculations, which is actually eliminated by including 
higher-order terms, m always remains a meaningful char- 
acteristic energy scale for strong coupling (see Appendix 
A). This is significant of a Mott insulating state, with a 
finite charge gap equal to 2m. 

An important point that can be extracted from the 
RG-flow is that the explicit divergence of the Umklapp 
scattering g u at the opening of the Mott gap results for- 
mally in: 



The jump of the LL exponent at finite T is significant of a 
Kosterlitz-Thouless transition, that is also accompanied 
by a sharp jump in the charge compressibility and Drude 
weight ||. Likewise, we like to stress on the fact that 
the strong decrease of the LL exponent by Umklapps 
has an enormous consequence on the small single-particle 
interchain hopping as well. 

Consulting Eqs.(]9|),(pT|), indeed we observe that this 
should also produce the complete vanishing of tj_ exactly 
at the Mott scale (and not only a partial reduction p5| ) : 



tj_(l > l e )=0. 



(12) 



In agreement with Giamarchi's assertion, we obtain that 
the Mott gap renders the single-particle hopping com- 
pletely irrelevant Q . We like to interpret the total van- 
ishing of t± at finite temperature as a sign of strong con- 
finement (of coherence) along the chains. For a physical 
explanation, do see next subsection. The divergence of 
intrachain Umklapp scattering produces in that case an 
insulating state with two degenerate bands. 

Now, we like to repeat that in absence of Umklapp 
scattering i.e. away from half-filling, this condition of 
confinement is never satisfied for the 2-leg Hubbard lad- 
der 1 14 1 . Precisely, to make interchain hopping irrelevant 
away from half-filling, one should make its amplitude de- 
creasing upon renormalization. This must be determined 
from inequality d± > 2 with 
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K> 



(13) 



K{1 > l c ) = 0. 



(11) 



the scaling dimension of the i^-perturbation. To obey 
this standard criterion of irrelevance || , one should start 
with a very small bare LL exponent K < 3 — V8 ~ 
0.171...: It is worth to note that this is never realized 
for purely on-site interactions and U <C t |l{|. This al- 
ways produces deconfinement away from half-filling that 
is characterized by a finite Zeeman-like splitting of degen- 
erate bands. The magnitude of the splitting A is defined 
in Eq. (|l|) . The system behaves as a Luther- Emery liquid 
p6| with dominant d-wave pairing. For a brief descrip- 
tion of the corresponding phase, see Section V. 

Nonetheless, the confinement condition away from 
half-filling might be performed in presence of long-range 
interactions where the chains behave rather as Wigner 
crystals fl5| . The underlying model in each chain is still 
of LL type, but the system develops a real tendency to a 
periodic arrangement or a 4fcp CDW: This hinders con- 
siderably the single-particle interchain transfer. It should 
be remarked that in that case, confinement might take 
place only at zero temperature. 

We also like to insist on the following point. Starting 
with very small interchain hoppings (such as m » fi) 
ensures g* = at zero temperature: Indeed, for the ID 
Hubbard chain and repulsive interactions it is well-known 
that the spin backscattering is (marginally) irrelevant 
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and then it can be omitted || . Note already the funda- 
mental difference with the deconfined regime, where spin 
backscattering also flows to strong couplings at the D- 
Mott transition producing immediately spin-gapped ex- 
citations 

B. Breakup of electrons on the chains 

The impossibility to get a finite single-particle inter- 
chain transfer at absolute zero might be naturally in- 
terpreted as follows. The system suppresses the bare- 
electron interchain hopping due to the strong constraint 
to avoid a double on-site occupancy, that is here explic- 
itly induced by the presence of the finite Mott gap. Now, 
we really like to emphasize that the vanishing of ij_ im- 
mediately at the Mott transition "caches" more a marked 
spin-charge separation or the complete breakup of phys- 
ical electrons. 

When T — > m, one already finds that the electronic 
Green function (e.g for Right movers) 



\J (vft — ix){uT — ix) 



(«r) 2 ]-" /2 , (14) 



tends to vanish for quite short time because the anoma- 
lous dimension of the electron increases fastly: r\ = 
l/4(K + K- 1 )-l/2 > 2 through the factor R- 1 ^), As 
usual, r denotes the Matsubara Imaginary time. Now, 
we show precisely that below the Mott transition there is 
no way to recombine a physical electron along the chains: 
This gives a simple explanation on the total vanishing of 
the bare-electron interchain transfer for T < m. 

First, from Eqs.©, one can observe that the velocity 
of charge excitations drastically decreases to zero. This 
tends to push forward the idea that charge and spin de- 
grees of freedom get really independent. Precisely, the 
opening of the Mott gap transforms charge excitations on 
the chains into fermionic Kinks, corresponding to pairs of 
doubly and empty occupied sites. The spin sector is still 
described by spinon- (anti)spinon pairs of the Heisenberg 
model. Do consult Appendix A for more explanations 
and details. 

Q Q 

In the following, the symbols T^" and S v * refer to a 
fermionic Kink and a (single) spinon, respectively. These 
are given by: 



Fpf » expiVnQt(-P®jc + Qjc), 



(15) 



Note that = ± refer to "electron" and "hole" like 
excitations precisely, and the spin of the spinon obeys 
S z = Qf /2 = ±1/2. The spinon and charge objects have 
different wave-vectors; an empty or doubly occupied site 
has a (double and opposite) wave- vector (—pQf)2kp. 



Second using Eq.(^), it is advantageous to decompose 
the bare electron operator 



rft - \r+ <?± 



(16) 



where describes the holon of the Fermi gas psfl : 



C f =eW£Q c ± ( 



-p$ ic + e jc ). 



(17) 



Now, it is sufficient to note that starting with a bare 
LL exponent K — > 1 (i.e weak U), one can equate: 



r + — \T + 



(18) 



This may simply explain why the recombination of elec- 
trons is definitely forbidden at the Mott transition: In- 
deed, this would require a fractional number of fermionic 
Kinks. As a natural consequence, the single-particle in- 
terchain transfer is already destroyed by the enhanced 
spin-charge separation arising at the Mott transition. 
This constitutes the main difference with the confine- 
ment phenomenon in Wigner crystals that occurs only 
at absolute zero. 



C. Prevalent fluctuations at the Mott transition 

Using known results on the Hubbard chain at half- 
filling and those of Appendix A, we get that the 4k f 
CDW and spin-spin fluctuations at q = it are then promi- 
nent around the Mott scale with the most diverging sus- 
ceptibility. 

Here, these are described by the operator: 



0^=^-+H.c, 



(19) 



that yields a non-zero expectation value below the Mott 
transition, resulting explicitly in: 

< Oi d Jx)Oi dw (0) > = < Oi d Jx) > = const., (20) 

and by the staggered magnetization operator: 

m i =<S^. £ r(<S^.)*+H.c. (21) 

The spinon operator has the quantum scaling dimension 
1/2, that produces: 



< mj (x) m3 (o)> = (-ir/x. 



(22) 



Since we have a jump in the LL exponent at the Mott 
transition, we recover Heisenberg correlation functions in 
each chain. 

Again unlike in the deconfined picture &0, here the 
opening of the Mott gap does not produce any spin gap. 
The single chain behaves as a pure Heisenberg spin chain 
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at the Mott transition. Our picture below the Mott en- 
ergy, however, is still incomplete. The tendency towards 
suppression of single-particle transport in the transverse 
direction is not the only effect of t± |^2|. The Hamil- 
tonian must be supplemented by extra "relevant" terms 
which are inevitably generated in the course of renormal- 
ization, by expanding the partition function as a function 
of t±_ ]l8[ |. For a detailed review on this important point, 
do see Ref. [||. 

IV. TUNNELING PROCESS 

It is maybe appropriate to rewrite explicitly the bare 
forward and backward hoppings 1 18 |2^] : 



7Y_l = t± cos v^Oe cos V^Os [ cos \A"^ C cos V^^s (23) 
+ (-l) x cos \fir§+ cos \/7r$+] , 

combining the boson fields in the two chains into a sym- 
metric "+" and antisymmetric "-" part. On the other 
hand, using Eqs.(Q) and ( [l6| ) one can also write: 



H± = £ t±e l{p '- p)kFX [c+iC;, 2 + H.c.] 
p,p' 

x [s+S-, 2 + R.c]. 



(24) 



Indeed, one always may interpret the hopping of a bare 
electron from one chain to the other as the hopping of 
a pure semionic charge = +1 and a single spinon 
with spin S* = ±1/2: The two objects behave then as 
half-electrons. 

As a consequence, we get the important identifications: 



CiC 2 + H.C. = COS y / TT& c COS ) 



f c~ 

and similarly for the spinon part 



(25) 



C^C_^ + H.c. = cos Vtt& c cos 1/7?$+, 



A. Spinon pair hopping 

As shown before, the total vanishing of t± for T <C to 
must be closely related to the breakup of electrons on 
the chains. On the other hand, we insist on the fact that 
"some" tunneling processes in t± 2 - respecting spinon 
and charge fermionic Kink separation - are still allowed 
in the far Infra-Rcd for T < to. Here, we like to empha- 
size that in previous works spin excitations at very low 
temperature have not been analyzed in great details p8| , 
or relevant tunneling processes have been completely un- 
dervalued |§. 

Now, we derive these properly as follows. Note that 
below we implicitly drop out oscillatory terms that do 
not influence the physics at the fixed point. 



First, it is important to mention that at the Mott tran- 
sition (I — » l c ), the relevant contributions are essentially 
furnished by: 

(C+C; 2 + H.c.) (c+ x C-„ 2 + H.c.) = cos (26) 

(C+CZ P2 + H.c.) (c+ +H.c.) = 00SV&+ 

Such operators acquire indeed a non-zero expectation 
value (see Appendix A): 



< cos \fAm<§~ > = < cos V^rr^t > 
whereas one still gets: 

< cos V47rO^ > = 0. 



TO, 



(27) 



(28) 



These may be primarily interpreted as tunneling pro- 
cesses of the massive fermionic Kinks lying along the 
chains. For a bare LL exponent K — > 1 one precisely 
gets: 

+ H - c -] 2 = [(^2) + H.c.] 271 , (29) 

= to[(C+C p - 2 )+H.c.] 2 . 

Approaching the Mott transition, second it is important 
to note that the (finite) amplitude of the relevant pair 
hopping(s) is given approximately by [See Appendix A 
and Eq. (A20)]: 

[t±(l -> lc)?/E F 2 « t 1 _ 2 exp2l c /E F 2 = i_L 2 m- 2 . (30) 

Ep is the Fermi energy. Now, averaging on the charge 
sector for I w l c , one gets the extra relevant Hamiltonian: 



SH in t = t^/mlY.S^S-^ + H.c.]' 



(31) 



p,p' 



SHint produces processes of coherent interchain spinon- 
(anti)spinon hopping triggered by the single particle one. 

Using the classification scheme of particle-particle and 
particle- hole hoppings from Ref. jig], one can equiva- 
lently rewrite M : 



SW-int = 35 cos V4tt$ s + (51 + .94) cos V 47re s (32) 
+ as cos \/47r$+. 



where we have: gi = t± 2 /m. Moreover, from Ap- 
pendix A, it is then clear that this generates an interchain 
Heisenberg interaction: 



SHint = J_l mi • m 2 , 



Jl = t ± 2 /m < 



(33) 



Since the spinon spectrum (driven by the kinetic part) 
is known to be invariant tuning the Hubbard interac- 
tion from weak to strong interactions, therefore the low- 
energy model is equivalent to two weakly- coupled Heisen- 
berg chains |20|. The intrachain Heisenberg coupling is 
here equal to the Fermi velocity, vp. 
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This proves rigorously that the interchain hopping is 
well sufficient to generate the antiferromagnetic spin ex- 
change J± « t± 2 /m > 0, that will make all spin excita- 
tions gapful. In particular, semi-classical considerations 
allow to predict the pinning of the 6 J and <&+ spin-fields 
at the fixed point (in the minima of the cosine potentials). 
Excitations will be separated from the ground state by 
an energy gap ~ J±. This leads to a ground state that 
would not have still a strong resemblance to that of the 
half-filled ID chain. Magnetic excitations and doping ef- 
fects are indeed very different. 

B. Effective Heisenberg ladder 

Excitations of the final spin Hamiltonian (H* +5Hi nt ) 
have been studied in great details in Refs. HtO|]. Rather 
than belabor the derivation, we only indicate the results. 
Rewriting the model in terms of four (real) Majorana 
fermions £ = 5Zi=i & an( i Pi one § e ^ s the following con- 
clusions. 

As for a general S=l magnet, the triplet excitations in 
a rung namely £ will have a gap - the Haldane gap |52| - 
m t = J±. The singlet branch p is located at quite high 
energy, of the order of m s = — 3Jj_ pi]] . For T <C J±, 
the asymptotic correlation function for spins on the same 
chain is given mainly by 

< m^m^O) > = (-l) x x -1/2 exp(-mta:). (34) 

This definitely produces a short-range Resonating Va- 
lence Bond solid exactly as in the deconfined limit. 

Note that taking only into account the (bare) inter- 
chain hopping term at a short wave vector q — 0, some 
spin excitations would remain gapless at the Infra-Red 
fixed point: 

SHint — 95 cos VAn<&~ + gi cos v4nQj ■ (35) 

The (s — ) modes are protected by the duality symme- 
try under $~ 8j, the (s~) sector would be in fact 
a critical point of the two-dimensional Ising type. The 
resulting Hamiltonian must be assigned to the same uni- 
versality class of the purely forward scattering model 
considered by Finkel'stein and Larkin Q and later by 
Schulz [Q , resulting in gapless- spinons and singlet (Ma- 
jorana) fermions. These excitations can be recombined 
and rewritten as a massless triplet magnon excitation £ 
in a rung. 

But, in the two-chain problem away from half-filling, 
again it is known that for general repulsive interactions, 
there is always a spin excitation gap and d-type pairing - 
or exceptionally orbital antiferromagnetism fluctuations. 
Here, the presence of such enigmatic spin liquid with rem- 
nant gapless magnon modes is definitely forbidden by the 
bare backward hopping at large wave vector q = tt. 



To summarize, we really want to emphasize that unlike 
in the deconfined region, the latter plays a central role in 
presence of strong confinement along the chains. 

C. Crossover to the D-Mott state 

We can now conclude that the two weakly-coupled 
Hubbard chain model is in the same phase C0S0 for very 
small bare interchain hoppings (confined region) and for 
rather large bare interchain hoppings (deconfined phase); 
CnSv denotes a state with n massless charge and v mass- 
less spin modes. This naturally demonstrates that the 
confincmcnt/dcconfincment transition at absolute zero 
in the half-filled two-chain problem is a simple crossover 
Now, let us compare symmetries of excited states 
in these two regimes. 

In the confined regime, we have a total spin-charge sep- 
aration on the chains. Charge and spin excitations are 
gapped but ruled by different energy scales. Addition- 
ally, we have shown that far below the Mott transition 
excitations can be classified as charge fermionic Kinks 
lying mainly along the chains with symmetry U(l) and 
charge ±1, and magnon like excitations with an under- 
lying symmetry SU(2) 2 x Z 2 @. 

In the deconfined regime, the effective model yields 
rather an enormous global S0(8) symmetry that can be 
briefly understood in the band basis, as follows P,pT|. 
The low-energy physics depends on a single effective cou- 
pling constant g (of the order of U) . In the band picture, 
the interaction part takes the specific form: 

~H ln t = -gcos\/47r$ Q cos\/47r$ ;) , (36) 

with (a, b) = 1,2,3,4 and bosonic operators are given by: 

($,6)! = ($,e)+ (*,e) 3 = (*,e)+, (37) 
(*,e) a = ($,e)- ($,e) 4 = (e,$)-. 

The indices ^ and ^ refer to symmetric/ antisymmetric 
charge and spin fluctuations in the band basis. It is 
then appropriate to use the re-Fermionization procedure 
(again, Q~ = -): 

* pa w expiVrrQ" (-p$ + 6 a ). (38) 

Then, for weak U the result is 

4 2 

H = Ys*Uv F T z d x t> a -g(*iTyy a ) . (39) 

a=l 

Pauli matrices r act on right and left sectors and ^> a = 
(*+o, *-a)- This is known as the S0(8) Gross-Neveu 
model. The latter has a remarkable property of "tri- 
ality" that is useful to equate various excited states in 
the deconfined phase. A remarkable fact that can be 
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shown from integrability of the model is that lowest ex- 
cited states are magnon like excitations with spin-1, and 
spinless charge ±2 or Cooperon and hole-pair. They 
are associated with two Dirac fermion excitations. This 
beautifully demonstrates preformed pairing with an ap- 
proximate d-wave symmetry (precisely, with a relative 
sign change between bonding- and antibonding pairs) in 
the deconfined state. Unlike in the confined region, spin- 
charge separation is here violated. 

Finally, note that in the crossover region i.e. when 
A w to the spin gap J± tends to rejoin the charge gap m. 
Both become of the same order as M oc exp — itvf/8U, 
the unique energy scale in the delocalized D-Mott phase. 
From the above analysis, one can furthermore predict 
that thermally activated fermionic Kinks on the chains 
with charge ±1 would then turn into excited Cooper pairs 
(in a rung). 



V. DOPING EFFECTS 

Doping the D-Mott state, it is well-known that the low 
energy excitations are only in the charge sector and cor- 
respond to a sound mode of the "hole pairs" . The phase 
is C1SO, i.e. a Luther-Emery liquid (26). The resulting 
system behaves then as a fluid of hard-core bosons with 
quite long-range d-wave pairing Pj7|Jll|] . 

Very close to half-filling, note that the SO (8) symmetry 
still remains [ p~T| , but this gets easily broken down away 
from half-filling in 50(6) x U(l). The U(l) symmetry 
comes from the hole-pair bosonic field that becomes 
critical, whereas the SO (6) symmetry results from the 
residual massive part. This is explicitly composed of the 
3 massive Dirac fermions ^2, ^3 and ^4 that are still 
coupled via a Gross-Neveu interaction. The most direct 
consequence is that spin-charge separation is only par- 
tially restored: There is a sixfold degenerate extended n- 
mode in the d-wave superconducting state that contains 
both massive spin and charge excitations Jli| . Note that 
unlike for the doped Hubbard chain, charge fluctuations 
in the superconducting state can only subsist at 4fcp [ |35| . 

It is important to remind that away from half-filling i.e 
in absence of Umklapp scattering, the 2-leg Hubbard lad- 
der is therefore described by a C1SO phase with promi- 
nent d-wave pairing and unconventional 5*0(6) x £7(1) 
symmetry. 

It is noteworthy that the hole-pairing phenomenon 
should subsist in strong magnetic field; The resulting 2- 
band model (of spinless fermions) predicts indeed promi- 
nent p-wave superconductivity for quite large t± J36[|. 



A. Hole doping effect on the confined state 

For completeness, now we study light doping effect on 
the confined state properties. 



Qd = \j- I dx d x ($ lc + $ 2c ) 



We remind that the Umklapp interaction g u has still 
to be taken into account for low dopings, S < to. It is 
therefore advantageous to take the half-filled picture that 
has been derived before, as a natural starting point. 

As usual, we model the doping by adding a chemical 
potential — /iQd with the charge operator: 



(40) 



2 y — / dx d *®c- 

Each chain is supposed to be equally doped. Note that 
the spin spectrum should not be affected by the low dop- 
ing effect producing inevitably the pinning of the fields 
and 0~. The resulting model then reads: 

H c=^J dx {(d x $+) 2 + (d x St) 2 (41) 

/K T 
— <9 X <I+ \. 

This Hamiltonian, describing a (purely ID-quantum) 
commensurate-incommensurate transition, has been ac- 
tively studied in the litterature. For a review, consult p. 
172 of Ref. ||. Formally, solutions of the equations of 
motion can be still written as fermionic solitons: 



(42) 



The bottom of this band is typically at energy-scales close 
to /i c = —to. Using Jordan- Wigner duality, one can also 
associate the hard-core boson field ||: 



A ps exp — iv^Oc 



(43) 



Using the fact that the fields and Oj are still pinned, 
then one can easily check that such object describes well 
a hole-pair with charge Qd = — 2 and zero momentum 
(Consult Ref. g p. 280): 



(44) 



Since the antisymmetric charge sector is not affected by 
the chemical potential, fluctuations of charge density in 
each chain become now strongly correlated: 



< d x <Pi c > = < d x <S> 2c >, 



(45) 



inducing pairing between holes of the two chains. Phys- 
ical excitations are rather described by the pairing field 
A that carries zero momentum (instead of T~ + ) . We like 
to emphasize that in this case, charge objects at- and 
close to half-filling are then different (fermionic Kinks 
with Q^r — ± turn into hole-pairs). 

Furthermore, one gets the important equality: 



[Q d ,A] = -2VKA. 



(46) 
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On the other hand, the (vertex) chiral operators J- p+ or 
A must also satisfy fl2ql: 



£„. = expi 



[Qd, A] = Q d A = -2A. 



Therefore, at the commensurate-incommensurate tran- 
sition one checks that the LL parameter K — * 1, ex- 
actly like for the lightly doped D-Mott state @. This 
reinforces the idea that at filling, there is no real differ- 
ence between the confined- and deconfined phases (T=0): 
They react similarly by doping. Such a universal value of 
the LL exponent in the two-chain problem has been first 
conjectured in Ref. and recently reached numerically 
in the strong U-limit with Density Matrix Renormaliza- 
tion Group approach J3£| (where the confined picture can 
be extended) . The limit K « 1 is consistent with the pic- 
ture of a very dilute hole-pair gas. More generally, this 
describes well a hard-core boson gas in low-density limit 
P0| . Note the difference with the one-chain case . 

One can also understand the preceding result as fol- 
lows. The holes stay in the same rung because they 
do not break spin singlets. The total charge mode 
becomes massless jl8 resulting in an expected Luther- 
Emery liquid (C1S0) [2§. Similarly to the lightly doped- 
Heisenberg ladder g2g or D-Mott state one im- 

mediately recovers prominent d-wave superconductivity 
with the pairing correlation functions: 



< A+(x)A(0) > 



-1/2 



(48) 



Note that adding eventually a tiny Coulomb repulsion 
interaction u between the chains, the triplet mass evolves 
slightly as m t = Jj_ — u, but remains finite H. 

Finally, for larger doping i.e when 8 > m, the intra- 
chain Umklapp channel can be neglected resulting in to- 
tal deconfinement, and the phase C1SO acquires an en- 
larged £"0(6) x U(l) symmetry and a sixfold degenerate 
"7r-mode" as well. Note that the hole-pair mode can be 
equally described in terms of the ©J or 0+ phases. 



B. Incoherence of t± away from haf-filling 

Finally, we would like to discuss the fact that the Lut- 
tinger model is not a good fixed point of the 2-leg Hub- 
bard ladder away from half-filling. This can be inter- 
preted as an incoherent effect of t± because the conse- 
quent renormalization of the interchain hopping away 
from half-filling turns the weak-coupling LL behavior 
onto a strong-coupling Luther-Emery fixed point. We 
may understand the incoherence of t± mathematically, 
as follows. An alternative definition of incoherence of 
the single-particle interchain transfer in LL's has been 
given in Ref. Q. 

For instance, one can check that charge eigenstates of 
the Luttinger Hamiltonian [^8| (Consult Appendix B) 



_i£ 
K 



e. 



(49) 



(47) are not (exactly) equal to Cr:-. As a result, the consequent 



renormalization of t±_ progressively hinders the coherence 
of anyonic-type excitations lying along the chains. The 
resulting competition between the interchain hopping tj_ 
and the (bare) interaction U produces inevitably a strong 
coupling fixed point with no charge object with charge 
Qf = ±1 (but paired holes) and then no spinon (but 
magnons). 

More generally, for doped N-leg Hubbard ladders (N > 
2) this also results in a complicated strong coupling fixed 
point with new stable excitations [Q. A LL fixed point, 
however, can still arise in some specific models e.g in 
the 3-leg ladder (or in the N-leg ladder, with N odd) 
very close to half-filling [[t2| or for a very anisotropic and 
restrictive network built with an infinite number of "very 
strongly" coupled chains JljJ. 

Finally, t± gets coherent only in the non-interacting 
case {K — 1). Spin and charge excitations have the same 
velocity and then eigenstates in each chain are usual Lan- 
dau particles: 



jf \ r + oil 

u pja — l'~P3 C 'pj\K=V 



(50) 



The ground state is simply achieved by making symmet- 
ric and antisymmetric combinations of the particle op- 
erators in the two gas. The fixed point is still a Fermi 
gas. 



VI. CONCLUSION AND DISCUSSION 

To summarize briefly, we have shown that the weakly- 
coupled two Hubbard chain model is in the same phase 
COSO, for irrelevant interchain hoppings (confined region 
with SU (2)2 x Z2 x [C/(l)] 2 symmetry) and for relevant in- 
terchain hoppings (deconfined region with S0(8) symme- 
try). At absolute zero, the confinement/deconfinement 
transition - again, occurring when the single-chain Mott 
gap is of the order of the Zeeman-like band splitting en- 
ergy - is a simple crossover plj . This statement i.e the 
absence of a metallic phase at half-filling, can be gener- 
alized for N-leg Hubbard ladders with N=3,4... (N is not 
too large). 

First, the phase of perfect confinement is always equiv- 
alent to a N-leg Heisenberg ladder. Taking into account 
spinon-pair tunneling process between successive chains, 
it is indeed sufficient to induce a pure spin-ladder with 
open boundaries in the transverse direction. Based on 
known results on coupled spin chains, we predict an in- 
sulating COSO phase for N even, and C0S1 for N odd [§. 

Second, similar conclusions can be reached for the N- 
leg ladder far in the deconfined regime (t± large), ana- 
lyzing the resulting N-band model fl^|. Let us briefly 
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insist on the main arguments. At half-filling, one gets 
the following Fermi velocities: 



Vj = v 5 = 2t{l - (^)"cos 2 (^ T )} i/Z , (51) 
with j — N + 1 — j. This results in: 



nj N I 1/2 



VN < V 2 = VN-1 < 



(52) 



Integrating properly the numerous RG-equations, one 
finds then a decoupling into band pairs and a hi- 

erarchy of energy scales 



Ma 



exp —irvj/U, 



(53) 



where the band pairs scale successively towards the D- 
Mott state of the two-leg Hubbard ladder. For N even, 
all excitations are then gapped (the phase is COSO). For 
N odd, the remaining band behaves like a single chain 
at commensurate filling resulting in a C0S1 phase. As 
long as N is not too large, then the deconfined region 
is still insulating and there is no fundamental difference 
with the confined part, where in contrast t± is strongly 
suppressed: The deconfinement/confinement transition 
is still a crossover. No metallic phase arises for large t± 
in the non-doped case. 

This completely leaves open how the crossover should 
evolve for an an increasing number of coupled chains 
to give back the metal/insulating transition observed in 
Bechgaard materials (N is very large). In particular, 
the decoupling of band pairs breaks down for large N 
(vi V2~-) and the analysis of the RG-fiow becomes very 
subtle mainly due to the plethora of relevant interband 
coupling channels Jl2[ . The fixed point for quite large 
t± is not yet known. A simplified (and much less rigor- 
ous) route to describe the strange metallic state of the 
TMTSF salt would be to start from the confined picture 
(with explicit spin-charge separation) and then to inter- 
pret the relevant hopping as an induced "self-doping" 
on the chains, that are not obligatory equally doped ||. 
Self-doping would appear as a possible way to mimic the 
small deviation of the commensurate filling due to the 
warping of the Fermi surface perpendicular to chain di- 
rection. 

Since here chains are not equally self-doped, then one 
does not predict any hole-pairing effect: Charge excita- 
tions would be then those of the single-chain problem i.e. 
doubly and empty occupied sites. For small induced self- 
doping i.e. for tj_ close to m, the resulting system might 
have similitudes with the lightly doped Hubbard chain 
This simple picture could explain the very small 
spectral weight 8(w) (Most of spectral weight comes from 
charge-gapped excitations) and the unusual frequency- 
dependent conductivity in the finite frequency regime 
p5|, observed on the basis of optical measurements |10|. 



On the other hand, it is important to remind that the 
TMTSF-family yields a good T 2 -resistivity showing the 



relatively large importance of the transverse hopping in 
this system. The induced self-doping would not be so 
light in these materials and then the ID LL-prescription 
based on a single doped Hubbard chain naturally breaks 
down (We remind that LL behavior is rarely stable in 
the N-chain model away from half- filling |44[). A possi- 
ble way to reconcile Fermi-liquid behavior with optical 
datas would be to study the effective infinite-band model 
where bands with intermediate indices would be suffi- 
ciently doped whereas bands with small or large indices 
would remain insulating. As soon as bands with interme- 
diate indices are consequently doped, one indeed expects 
a weak-coupling fixed point of Fermi-liquid type |ff4| . 



ACKNOWLEDGMENTS 

We acknowledge discussions with L. Degiorgi, F.D.M. 
Haldane, U. Ledermann, A. Nersesyan, R. Noack, T.M. 
Rice and thank T. Giamarchi for a correspondence. 



APPENDIX A: EXCITATIONS AT THE MOTT 
TRANSITION 

Note that the opening of the Mott gap produces charge 
fermionic kinks confined to the chains. 

Using the transformations (K « 1 denotes the bare LL 
exponent) 



'I',, = \ l\'\>,. & jc = (-).,. \ K, 



(Al) 



the interacting part driven by umklapps in the chain j is 
of Sine-Gordon type: 



Hi 



-g u cos V 4irK$j C — —m cos V4n&j c . (A2) 



Solutions of the underlying Sine-Gordon model arc 
known as soliton-like; Solitons are fermionic particles: 



fpf = K P j exp i^Qf(-p$j c + <dj C ) 



(A3) 



with the chosen gauge K + jK^j = i, and Qf — ±1 re- 
fer to "electron" and "hole" like excitations respectively. 
Therefore, it is not surprising that the charge part of the 
Hamiltonian density can be refermionized, as [|| 

n c = J2 • '>'■'■):/>, :'■)„ - , + h.c. (A4) 



Excitations can be viewed as pairs of doubly/empty occu- 
pied sites; a Kink carries momentum {—pQ^yikp. 
It should be noticed that fermionic (exchange) statistics 
means precisely: 

^pj( x )^Pj(y) = F$j{y)7tj(%) ex P -hsgn(x ~ y), (A5) 
with 
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7 



(A6) 



Using Fourier and Bogoliubov transformations, the band 
structure of these solitons reads 



E(k) = ±y/{uk) 2 +TO 2 



(A7) 



This results in a semi-conducting picture, where the 
charge gap is equal to (2m). At zero energy, there is 
no way to create pairs of doubly/empty occupied sites. 
The ground state yields a long-range Akp CDW. 
To study pair-hopping processes, one can rewrite: 



(A8) 



Tiint = -2g u cos V47T$+ cos V47r$ 



At the Mott transition, the pinning of these terms renders 
one-point correlation functions constants 

< cos V&r&t > = < cos V47r$~ > Kim. (A9) 

We use normal ordering of phase exponentials in Sine- 
Gordon models: 

< cos/?* > = m K02/47T : cos /3$ :, (A10) 

and the fact that the pinning phenomenon produces: 



cos v 47r$ 



cos V 47T$ 



1. 



(All) 



Due to cluster decomposition principle, the associated 
two-point correlation functions become also constants at 
large distance. Note that the one-point correlation func- 
tions of the dual fields are still equal to zero (the 2-point 
correlation functions decay exponentially). 

At the Mott energy, spin excitations or so-called 
spinons tend to be also localized on the chains. This car- 
ries the spin of a physical electron and then is described 
by the operator: 



(A12) 



with, S z — Qf /2 = ±1/2. Such objects are known to be 
located at wave vectors pkp, and obey a semionic statis- 
tics with 7 = 7r/2 Pq| . In absence of charge fluctuations 
and for vanishing tj_, spin fluctuations in each chain are 
produced only by spinon pairs p2| . 

As for the antiferromagnetic Heisenberg chain, it is 
appropriate to define the staggered spin operator as: 

m, = S^a{S^)* + H.c, (A13) 
resulting explicitly in: 

rrij (x ( cos v / 2lr8j S , — sin v^Ojs, cos V^r'I'js). (A14) 

Such definition of nx, respects the total breakup of phys- 
ical electrons below the Mott scale. Now, it is maybe 



important to compute explicitly the product mi • m.2 
(Consult subsection IV- A). One gets: 



m i m 2 



2 cos\/27r(ei s - 9 2s ) 
-I- 2 cos V2tt$ 1s cos \/2tt§ 2 
= 2 cos \/47r07 -I- cos VItt<S> l 

^From Abelian equalities: 

H.c. = cos \/ttQ~ cos y/n$>J , 
H.c. = cos cos \/7r$j", 
this also results in: 

1 2 



(A15) 



COS V 47T$ 



mi • m 2 = [J2 S p1 S p'2 + H - c -] 



(A16) 



(A17) 



This gives an explicit relationship between antiferromag- 
netic interchain coupling and spinon-pair hopping. 

Now, let us comment on the amplitude of the spinon 
pair hopping(s) in the entrance of the Mott transition 
i.e. I — > l c [See Eq. (31) in the Section IV]. At very short 
distances I — > 0, pair hoppings are neglectable. Under 
renormalization of the short-distance cut-off, the associ- 
ated amplitude is known to obey [fL8[: 



^ = (l-K)g + z(l/K-K), 



(A18) 



with the bare conditions for I — > 0: g(0) = 0, z(0) = 
t ± 2 /E F 2 and (1/K - K) ~ U/v F for small U. E F is 
typically the Fermi energy. See e.g. Ref. Q page 226. 
For / < l c , this can be simplified as: 

dl 

The amplitude of the pair-hopping (s) becomes finite ap- 
proaching the Mott transition. Using Eq. (0), we find: 



z(l) x U/vp- 



(A19) 



g(l c )oc z{l)dl= [ ° dz{l)fat ± 2 /m 2 
Jo Jo 



(A20) 



We have neglected z(Q) in front of z(l — * l c ) — 
t±(l — * l c ) /E F 2 = tj_ 2 /m 2 because m -C Ep. Further- 
more, averaging explicitly on the charge sector, we get 
that for T — * m the effective amplitude of the spinon 
pair-hopping(s) reads: 

9i{h) = §(lc) < cosa/47t$+ > cx t±_ 2 /m. (A21) 

Below the Mott transition (I > l c ), t±(l) and K(l) are 
both renormalized to zero [See Part III, Eqs. (11) and 
(12)]. Therefore, <?i now obeys: 



dgi 
dl 



9i, 



(A22) 



with the new bare condition gi{l c ) = t± 2 /m. Far in the 
Infra-Red, the spinon pair-hopping(s) will "diverge" , pro- 
ducing a spin gap of the order of J± = t± 2 /m in the 
spectrum (which can be exactly obtained via refermion- 
ization, see Part IV B). 
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APPENDIX B: CHARGE EXCITATIONS IN A LL 



In a LL, the charge Hamiltonian is plasmon-like: 

K = h\ dX lt {Pc ~ Po)2 + K( ^ e ^- ( B1 ) 

dx^c = (Pc ~ Po) measures fluctuations of charge den- 
sity, and VO c is the conjugate momentum to $ c . The 
equations of motion are usual d'Alembert equations. 
It is appropriate to use the chiral decomposition: 



K 



(B2) 



Again, p — ± refers to the direction of propagation (right 
or left). The associated Hamiltonians are defined as 



H ± = ^-Jdx (ve±) 2 



The objects 0± are chiral i.e. they obey p8[ 



[e ± (x),T^d v e±(y)] =iS(x-y). 



(B3) 



(B4) 



Chiral vertex operators of the effective Gaussian model 
read 



Cl a = expiJ^Qfe ± . 



(B5) 



They describe charge excitations of a LL. In particular, 
these have a charge Qf — ±(1) because 



[Qc, Cp " ] — Q c Cp " . 
The charge operator is correctly normalized as: 



Qc 



dx d x $ c . 



(B6) 



(B7) 



Such objects correspond to anyonic excitations since they 
obey anyonic commutation relations with: 
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2K 



(B8) 



Note that the Hubbard interaction between physical elec- 
trons produces (at low-energy) a change in the statistics 
of "holons". For the free electron gas, one gets K = 1 
and then holons are rather scmions or half-electrons. 
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